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We further analyze the recently proposed criterion, according to which a superfluid transition occurs 
whenever the r.m.s. displacement of a particle from its initial position in imaginary time reaches a 
certain fraction of the interparticle distance. The critical temperature T\ is expressed in terms of 
the single particle mobility. Within the Drude approximation for the mobility T\ is determined by 
the friction coefficient (or by the viscosity for dense liquids). To check the validity of the criterion 
the resulting formula is applied to non-ideal Bose gas, to liquid helium and to 3 He- 4 He mixtures. 
A reasonable qualitative agreement with available experiments is obtained in all cases. In the case 
of the Bose gas the present approach results also in an upper bound on the peak value of Tx/T®, 
where T° is the critical temperature for the ideal gas. 



PACS numbers: 67.40.-w, 05.30.Jp, 67.20.+k 
I. INTRODUCTION 

Recently considerable attention has been attracted to 
the long standing problem of how the critical tempera- 
ture of the superfluid transition depends on densjte, (or 
interaction strength) in interacting Bose systemsoili. If 
the interaction is switched off then the temperature of the 
Bose-Einstein condensation is given by the well known 
ideal gas formula 
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where n is the density and m is the particle mass (we use 
units Kb = 1). This formula should be approximately 
valid also for an interacting system in the limit of low 
density. By now it is generally believed that if the den- 
sity increases T\ first rise with respect to the ideal gas 
value, but the exact form of the correction to Eq. ([!]) still 
remains, to some extent, a point of contmversy, despite 
the large number of recent investigationsoEl 

While for a weakly interacting gas there exists (at least 
formally) a small parameter, proportional to the scatter- 
ing length, the problem looks even more difficult if the 
density is not small. The critical temperature should fi- 
nally goes down if the density is further increased, as ob- 
served e.g. in liquid helium, primarily due to an increase 
in effective mass, that is believed to replace m in Eq. (|l|) 
for an interacting system. Since T\ is obviously not an 
universal quantity (unlike critical exponents it depends 
strongly on microscopic details of the system) it seems 
reasonable simply to evaluate the transition temperature 
numerically for a given system, using-e.g. the Path Inte- 
gral Monte Carlo (PIMC) methodB'n, but some analyt- 
ical results, if available, would certainly lead to a better 
understanding. In case of dense systems with strong in- 
teraction a lattice model may provide a good starting 
point (see e.g. Refs. |^Jl4|) though the lattice introduces 
some additional commensuration effects not present in 
continuous systems. 

On the other hand a simple phenomenological super- 
fluidity criterion was proposed recently by the authorlla, 



which results in an analytical expression for T\ in a dense 
system. The approach is based on a Feynman picture of 
a bosonic system evolving in imaginary time, when the 
superfluidity is related to the appearance aLpermutation 
exchange cycles of arbitrary large lengthOEl At high 
temperatures only the identity permutation is important 
in the partition function and each particle returns to its 
initial position in the imaginary 'time' (3 = l/T . Its 
root mean square (r.m.s.) displacement R is essentially 
the corresponding quantum coherence length which in- 
creases when the temperature is lowered. The superfluid 
transition occurs when R finally reaches the value of the 
order of the interparticle distance a = nT 1 / 3, (at lower 
temperatures nontrivial permutations can be no longer 
ignored). This picture suggests that one can first evalu- 
ate R(T) at a given temperature T and then obtain an 
estimate for T\ e.g. from the condition R 2 (T\) = £a 2 , 
where £ is some numerical factor of order unity (see Sec. 
II for details). 

This approach is similar in spirit to the well known Lin- 
demann melting criterion, which states that a solid melts 
whenever the r.m.s. displacement of particles from their 
equilibrium positions (due to lattice vibrations) reaches 
a certain fraction of the nearest-neighbor distance. It is 
known that in many cases the ratio of these two quanti- 
ties (the Lindemann ratio) is almost universal thus mak- 
ing the criterion quite useful for qualitative estimates (see 
e.g. Ref. |l8|). One may hope that the present approach 
which treats the superfluid transition as some kind of 
'melting in imaginary time' may also be of use at least 
as a first approximation. 

From the physical point of view R may be regarded 
as a generalization of the thermal dc-Broigle wavelength 
At = {2itTi 2 /mT) 1 / 2 , which is the coherence length for 
non-interacting particles at a temperature T. Indeed, for 
the ideal gas R ~ At and our criterion coincides with 
the textbook argument which states that Bose-Einstein 
condensation occurs whenever At ~ a (which yields 
T\ ~ (h 2 /m)n 2 / 3 ). If the interaction is switched on, 
interparticle scattering will tend to reduce the coherence 
length, so that R can be much smaller than At- The 
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physics here is essentially the same as in the decoherence 
due to environmentEj, since for a given particle all other 
particles will now act as some dissipative environment. 
One may also recall that the kinetic energy K(T) of a 
particle in a quantum liquid is no longer given by |T 
but is higher due to interaction (at low temperatures the 
mean kinetic energy tends to a constant value determined 
by a zero-point motion), so that the particle wavelength 
estimated as ~ (h 2 /mK (T)) 1 / 2 indeed gets smaller than 
Xt- 

If the interaction leads to the reduction of the co- 
herence length it also reduces the superfluid transition 
temperature (the initial growth of T\ with interaction 
strength in a Bose gas has a different origirfl and will be 
discussed in Sec. III). Thus e.g. the observed decrease 
of T\ with pressure in helium may be understood either 
as a decoherence phenomenon or as a consequence of the 
more intensive zero-point motion in a pressurized system. 

In Ref. |lj this supcfluidity criterion was applied to 
molecular hydrogen H 2 and results in T\ ~ 1 K fa^thc 
temperature of the possible superfluid transitiorcjtj. 
This estimate is significantly lower than thep-very first 
one, based on the ideal gas formula (|]), 6^-8 KEj (though 
consistftat with lateij-estimatesEil and PIMC results for H2 
clustersE-3 and filmsca). This result was obtained, how- 
ever, by assuming that the 'Lindemann' parameter £ is 
universal, i.e. the same as in case of liquid helium. Such 
an assumption may cause some doubt and its further 
analysis is certainly in order. The r.m.s. displacement 
R that enters the superfluidity criterion was evaluated 
in Ref. |l3| using the Caldeira-Leggett mode]Eil which also 
needs some justification in our case. 

In this paper we study this problems in more detail in 
order to better understand the accuracy and reliability 
of the criterion proposed. In Sec. II we review some 
general aspects of the theory. The general formula will be 
derived, which relates R 2 to the mobility of an external 
particle in the liquid. We will show that the Caldeira- 
Leggett model approach is equivalent in our case to the 
Drude approximation for the mobility. Within the Drudc 
approximation the mobility is determined entirely by the 
collision rate 7 and it appears possible to derive a general 
formula for the superfluid transition temperature which 
is an universal function of hj/T®. If this ratio is much 
larger than unity the superfluidity is suppressed. 

Next we apply the formula obtained to different su- 
perfluid systems and compare the results to available ex- 
perimental data. In Sec. Ill a weakly interacting Bosc 
gas is considered. Here we analyze the ratio of T\/T® 
as a function of density and compare the result to the 
known PIMC calculations. Sec. IV is dpzpted to liquid 
helium. Contrary to the previous worktj where 7 was 
related to the mean kinetic energy of a helium atom, the 
transition temperature is expressed here as a function of 
the shear viscosity in the normal state (at T close to T\). 
We also check the validity of the Drude approximation 
and find out that for liquid helium this is quite accurate 
if exchange effects are neglected. Finally, in Sec. V we 



extend the present approach to 3 He- 4 He mixtures. With 
no additional fitting it appears possible to describe qual- 
itatively the behaviour of T\ as a function of both the 
total density and the concentration of 3 He atoms. 

II. GENERAL FORMALISM 

According to the Feynman's path integral approachE] 
the partition function Z for a system of N interacting 
bosons can be represented as a sum over all possible per- 
mutations of these particles 

z =ME/n*«/n Br '( T )«p(-J 5) ' (2) 

P J i J i 
rhf3 -2 

S= / [E^ + E^-r^r 

where V{ri —Vj) is the interparticle interaction potential 
and (3 = 1/T. The integration in Eq. (||) is over all 
paths with rj(0) = 17, Vi(hf3) — P17, where P is some 
permutation of ./V particles. 

At high temperatures only the identity permutation 
contributes to the sum, but as [3 increases, other per- 
mutations finally also become important and the system 
becomes superfluid. In this picture the superfluid transi- 
tion is related to the appearance of long cyclic permuta- 
tions (see e.g. Ref. |l^ for details). This means that the 
critical temperature should be roughly of the same order 
as a temperature at which nontrivial permutations can 
no longer be neglected. Then we may try to estimate T\ 
from a criterionEj 

R 2 ^— j ([r(r)-r(0)] 2 )dT = £a 2 , (3) 

where r(r) is a path of some arbitrary particle in the 
system of distinguishable particles and £ is a numerical 
factor, to be determined later. Average in this formula 
means the functional integration over all paths of the 
chosen particle r(r) with r(0) =r(/3) 

(...) = || Vv{r){...)^ V L\J % J dr^f-W{r(r)} 

(4) 

where the effective action W / {r(r)} arises from the inte- 
gration over other TV — 1 particles (i — 1,2 ... N — 1) with 

exp(-W/{r}) = J []pr l (r)exp(-^'{r,r l }) , (5) 
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The condition (g) merely states that near the transi- and hence 
tion the mean displacement (in imaginary time) of a given 
particle from its initial position becomes comparable to 
the interparticle spacing. This is similar in spirit to the 
well known Lindemann melting criterion. For this reason 
we shall some times refer to Eq. (||) as the 'Lindemann- 
like' criterion and also call £ the 'Lindemann ratio'. We 
shall assume here that £ is roughly a constant, but, cer- 
tainly, possible weak dependence of £ on system param- 
eters (e.g. on density) cannot be excluded and may be 
important in some cases (see e.g. Sec. III). 

In the ideal gas the left hand side of Eq. (|J) is essen- 
tially the square of the de Broglie thermal wavelength 
A T = 2irh 2 /mT, and our criterion looks like At ~ a, 



but interactions will tend to reduce R 2 (this was ob- 
served e.g. in Ref. ^). One may think of this reduction 
as arising from the decoherence due to interaction with 
environment. Neighboring particles which scatter from 
the one we are looking at 'measure', in a sense, its posi- 
tion thus reducing its quantum uncertainty in coordinate 
space (compare e.g. with Ref. |l9| ). 

Estimating R 2 in a system of interacting particles is 
still a very complicated problem, even if exchanges are 
neglected. One has first to calculate somehow the ef- 
fective action W^{r} given by Eq. (Q). However, since 
we are interested only in the simple average (g) we can 
avoid the evaluation of VF{r} and relate this average to 
the mobility of an external particle. 

If we introduce a Fourier transform of the path r(r) 
according to 



r(r) 



E 
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where w B = 2im/h/3, n — 0, ±1,±2, ... then in three 
dimensions the mean square displacement defined by Eq. 
mf) can be written as 
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Now let us write the expression for the mobility in the 
Matsubara representation /zm(w„). For this purpose we 
introduce an external force f(r) = f exp(— iui n r) acting 
on a given particle and a perturbation term 
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in the action. The force will induce the average veloc- 
ity (v(t)) = v„ exp(— iuj n T) and the mobility fiM^n) is 
defined according to 

V n =/iM(Wn)f- (8) 

Since (v(r)) = i(r(r)), to the first order in f we obtain 
(v(r)) = - jf (r(r)r(r'))f(r')dr', (9) 



A<m(w„) = uj„f3(\r ri 



(10) 



This means that the average we are interested in can be 
expressed as follows 



r 2 = — ^m(uju) 
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The Matsubara mobility is related to the mobility at real 
frequency n(co) by the usual equation 

Miw(^n) = Ki\uJn\)- 

Since fi(u>) is a response function and is analytical in 
the upper half-plane of the complex frequency one can 
easily verify that 



(i(iu n ) = - 

7T 
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for positive u> n . Substituting this expression in Eq. ( |ll| ) 
and carrying out summation over n we finally obtain 



— Re/x(w) 



w ?kt\ 2T 
coth( — ) 
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Thus if we know [m(u>) from some theory or experiment, 
we can evaluate the mean square displacement (|l3|) which 
enters in our superfluidity criterion. 

Important asymptotics can be obtained, however, 
without knowing the exact expression for n{u>). We shall 
assume only that the static mobility /1 = /i(0) is finite 
and introduce the scattering rate 7 = l/(m/x). First of 
all , if the interaction is switched off, then R 2 — > h 2 (3/2m. 
Then the criterion (g) results in T\ = (l/2£_)h 2 n 2 / 3 /m. 
Since this is the ideal gas formula for the critical temper- 
ature (|l|) we conclude that l/2£ ~ 3.31, i.e. 



£ ~ 0.15 



(14) 



If however the interaction is so strong that E7 > T 
then the leading low temperature contribution to R 2 will 
be logarithmic in frj/T 



i?2^ H^4p)+ const 



(15) 



It is clear from (|l3| ) that this logarithmic behavior is 
actually independent of the form of //(w) provided the 
static mobility is finite i.e. /i(0) 7^ 0, and the prefactor is 
determined entirely by /1 = 1/ (7717) . The second term in 
( |l5| ) depends on the exact form of the function (J>(ui). 

It is important, that though we consider low temper- 
atures, we always deal here with nondegenerate liquid 
where the mobility (and hence 7) is approximately inde- 
pendent of temperature (see Sec. Ill) . Below the A point 
the friction of an external particle is determined by col- 
lisions with quasiparticles (phonons and rotons) and will 
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tend to zero as T — > 0. For this reason it is more ade- 
quate to describe the system below the A point in terms 
of an effective mass m*. Our present approach is in fact 
complementary to the one based on the notion of effective 
mass, which results in T\ ~ (Ti 2 /m*)n 2 ^ 3 . The temper- 
ature dependence of R 2 in the normal system above T\ 
(with all exchange effects switched off) can hardly be de- 
scribed by a simple renormalization of the particle mass, 
and the incoherent diffusion may be more adequate a 
model for the particle motion (see also Sec. IV) . Though 
the effective mass approach results in the apparently sim- 
ple expression for T\ it is not really simpler, since there 
are no analytical closed expressions for m*. 

Using Eq. dD5| ) we can obtain from Eq. a very 
simple formulalij 
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where a is some unknown factor, dependent on the high 
frequency behavior of and To is essentially the crit- 
ical temperature for the ideal gas (it differs from Eq. (|l|) 
only by a numerical factor). 

Thus the temperature of the superfluid transition cru- 
cially depends on the ratio fvy/Tq which is a function 
of density and interaction strength. At hj <C To the 
ideal gas formula (Q) is valid while in the opposite limit 
the critical temperature of the A transition is exponen- 
tially small due to the decoherence phenomenon. Since 7 
should normally increase with density (for systems with 
strong repulsion between particles at small distances), 
the formula obtained qualitatively explains the suppres- 
sion of Tx m sufficiently dense systems. 

Now we introduce a simple model for At high 

frequencies (higher than the collision rate) the mobility 
is the same as for a free particle, i.e. ~ i/(mw). 

Then we may try to use, as a first approximation, the in- 
terpolation formula for the frequency dependent mobility 
in the Drude form 
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where A is the friction coefficient which determines the 
static mobility /1 = fi(0) = 1/A. Frequency dependence 
of this kind is valid for the mobility of a heavy particle, 
moving with friction according to the Langevin equation. 
Variational calculations of Ref. suggest that the Drude 
formula may be adequate for a light particle as well. In 
fact this interpolation formula is reasonable when there is 
only one frequency scale in a problem. Strictly speaking 
the friction coefficient should be also frequency depen- 
dent, i.e. A = \(u>), but in many cases it is possible to 
neglect this dispersion since Re/i(w) will be already very 
small when the dispersion in A becomes significant. Now, 
substituting Eq. (17) in Eq. (|13[) we obtain 
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where 7 = A/m, ip(x) is the digamma function and 
C ~ 0.577 ... is the Eulet's constant. This very expres- 
sion was used previously^ when the motion of a particle 
was described by the Caldeira-Leggett model with ohmic 
dissipation. Note that this formula (with the suitable 
choice of 7) may be used also to describe the suppression 
of superfluidity by an external disordered potential (this 
will be discussed elsewhere). 

The final result for the transition temperature in case 
of n(u>) given by Eq. (|l7|) may be written as 
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7 = A/m 



(19) 



where the function F(x) is implicitly defined by the equa- 
tion 



x = C + ip(l + l/F(x)). 



(20) 



From Eq. ( ^o| ) one can obtain F(x) ~ (6/ir 2 )l/x at 
x — > while F(x) becomes exponentially small F(x) ~ 
exp(C — x) at x — > 00. In the intermediate region x ~ 1 
this function is well approximated by the interpolation 
formula F(x) ~ 1.6exp(— 0.96x)/(l — exp(— 0.96a;)). 

Unfortunately, even for such a simple choice of fi(u>) 
there seems to be no opportunity to calculate the value 
of 7 analytically for a dense liquid with a strong in- 
teraction between particles. For this reason one has to 
use experimental results as well as approximations with 
some adjustable parameters. One can also try to use the 
fluctuation-dissipation theorem for the velocity fluctua- 
tions and write the following expression for the kinetic 
energy 



K(T) = ~mj ^ftu;Re M (u;)coth(^) 



If we put T = then 
3 

K = —m 



(21) 



(22) 



where K = K(0). This integral is known to diverge log- 
arithmically for the Drude mobility at high frequencies. 
If we assume that for true /i(w) there is only one charac- 
teristic frequency scale 7 then obviously 

K ~ mh<y 2 /\ ~ 717 

This approximation was used in the previous paperlll, 
and it has some advantage, because the kinetic energy 
K (essentially the zero-point energy) can be easily evalu- 
ated for various systemsEEI. In what follows, however, we 
shall try to make a more direct estimate of the damping 
parameter 7. 



III. WEAKLY INTERACTING BOSE GAS 

Before we turn to dense liquids consider first a non- 
ideal Bose gas. In the gas the transition temperature 
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is close to Eq. (||) and we are interested in corrections 
to T° at non-zero density. To estimate these corrections 
consider a classical gas of hard spheres of diameter d. 
Then the scattering rate is given by 



7 ~ v/l 



■nnd 2 v 



(23) 



where I ~ l/(nird 2 ) is the mean free path and v 
(T/m) 1 / 2 is the mean velocity. Since T 

K 2 n 2 / 3 /m 
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7 „ A d 2 n 4/3 



(24) 



If we introduce a small parameter z = {nd 3 ) 1 / 3 then we 
see that the ratio hj/T® ~ z 2 is small at low densities 
and one can expand the formula ( [L8| ) for i? 2 in fry/T. 
Retaining only the first correction in 7 we find from Eq. 
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The correction due to decoherence is negative and leads 
to the suppression of T\. It is generally accepted, how- 
ever, that, ptiic first correction to Eq. (Q) should be 
positiveta ~Ut3. The only way to account for this increase 
of T\ in the present approach is to assume that the 
parameter £ also slightly depends on density and gets 
smaller as the density is increased i.e. 



(26) 



where £0 — 0.15 and the small density dependent cor- 
rection 5 will be specified later. The physical reason for 
such a dependence is the same as in Ref. 0: particles in a 
more dense system tend to be more homogeneously dis- 
tributed through the whole volume, hence it is easier for 
a particle to find a neighbor for exchange. Then we can 
write 
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Now we need some more accurate estimate of 7. For 
a nondegenerate system of hard spheres one can evalu- 
ate the friction coefficient A = from the Bolzmann 
equation and in the low density limit 



A = ^nd 2 (TrmT)^ 2 



(28) 



(see e.g. Ref. g^). The additional factor of 4 is introduced 
here because in our case, contrary to the pure classical 
one, the wavelength of the particle is larger then d and 
the s-wave scattering cross section is known to be 4ird 2 
rather than nd 2 as for a classical sphere. 
Hence near the transition, when T 
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Most of recent calculations suggest, that at low density 
the first positiy*|-carrection in Eq. ( |27f ) should be of the 
first order in cfiHTa, so in this limit we can take 5 ~ z. 
At higher density the Lindemann patio is likely to reach 
some finite limit close to £ ~ 0.12.Lj Then we may try a 
simple interpolation formula 



5 = 



1 + a 2 z 



(30) 



Then the final expression for the transition tempera- 
ture can be written as 
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1 + a 2 z 



a 3 z 
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3U/3 



(31) 



where 03 ~ 1.2 from Eqs.(|27|) and (J29J), while ai and a 2 
should be taken from some microscopic theory or exper- 
iment. The dependence of T\/T® on density given by 
Eq. pi] ) is shown in Fig. 1 for gu — 2.3 as obtained 
e.g. in recent Monte Carlo studytl and for a 2 — 9.9. 
This value of a 2 was fixed by the requirement that at 
z ~ 0.61 (which approximately corresponds to helium 
density, when nd 3 ~ 0.23) £ should be close to CL12 i.e. 
S ~ 0.2. The squares represent the PIMC dataQ for a 
quantum system of hard spheres. Though we have taken 
rather large initial slope a\ ~ 2.3, which is almost an 
order of magnitude larger than reported in Rcf. 0, the 
resulting curve does not deviates much from the PIMC 
data (the largest deviation is less than ten percent) and 
reveals a qualitatively similar behaviour. 
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FIG. 1. Critical temperature T\ of an interacting Bose gas 
versus density from Eq. ( |3l| ) with a\ = 2.3, 02 = 9.9 and 
as — 1.2 (T° is the critical temperature of the ideal gas, d is 
the hapd-sphere diameter). For comparison, PIMC simulation 
resultsu are also included. 

What is perhaps more interesting, explaining the ini- 
tial rise of T\ with density by a change in the Lindemann 
ratio £ places an upper bound on the peak value of T\/T®. 
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Indeed, if the parameter £ changes from £o ~ 0.15 to 
something like £ ~ 0.12 at high densities then the maxi- 
mum value of S in Eq. ( p7| ) is 0.2 and hence T\/T° should 
not exceed 1.2. This result is consistent with the PIMC 
calculationsQ but is incompatible with recent experiment 
in 4 He-Vycor systemtJ where the peak value of T\/T® 
was found to be almost 3.5 times larger. This disagree- 
ment seems to be not yet explained and one may suspect 
that we do not fully understand all possible effects of a 
disordered porous substrate in Vycor samples on T\. 



IV. LIQUID HELIUM 

Let us see first what an elementary kinetic approach, 
can suggest for It is known from experimentE2l 

that above the superfluid transition the momentum dis- 
tribution of helium atoms is approximately of a Gaussian 
type with the mean kinetic energy K ~ 15 K. If we used 
a Boltzman equation for a classical gas of hard spheres of 
diameter d at effective temperature T* = |_K" then in the 
relaxation time approximation we obviously would obtain 
the frequency dependent mobility in the Drudc form (fiT 



A = Airnhn 



-1/3 



(34) 



with the inverse relaxation time 7 ~ v/l ~ irnd 2 v as m 



Eq. (|2^) where now 
one has 

A - 



(K/m) 1 ' 2 . At the same time 



7777 



Trnd 2 mv 



for the friction coefficient and 



77 ~ mvnl ~ (l/ir)mv/d 2 



(32) 



for shear viscosity. From these formulas we get 

A - Tr 2 (nd 3 )rid (33) 

Extrapolating this result to a dense system with nd 3 ~ 1 
we see that the friction coefficient in a liquid is in fact 
proportional to the viscosity. Note that the estimate (U4) 
is quite reasonable for liquid helium above the A pointES. 
Indeed, the measured viscosity above T\ is almost inde- 
pendent of temperature and is equal to 77 ~ 3.5 • 10~ 5 
poise (at s.v.p. and at T = 3.2 K) while the more ac- 
curate version of Eq. ([32]) 77 ~ 0.18(mT*) 1 / 2 /c? 2 gives 
77 ~ 3.6 • 10~ 5 poise at d ~ 2.2 A. 

Practically the same result as Eq. (|3^) arises also from 
the hydrodynamic approach. For a particle of a large 
radius 1Z 3> a in a liquid the friction coefficient A from 
the Stokes' formula is QnrfR,. It is clear, however, that for 
a particle of the liquid, which is not large as compared 
to other particles, this formula should be modified. Long 
ago it was argued (see e.g. Ref. [32] and references therein) 
that it this case one should rather take A = AirnlZ, as for 
an empty bubble in a liquid. One may suppose in this 
connection that a particle in a liquid is surrounded by a 
'correlation hole', so that one may view the particle as 
associated with a 'bubble' of radius 1Z ~ rT x l z . Then we 
have an estimate 



where we have put 1Z = bn~ x l 3 and b is some numerical 
factor of order unity. Certainly, this estimate makes sense 
only above the A transition and for sufficiently dense liq- 
uid. When nd 3 ~ 1, estimates ([}3|) and (|34| ) lead to 
essentially the same results. In what follows we shall use 
Eq. (|34|) since this will result in a slightly more trans- 
parent final formulas. 

Now, making use of Eq. ( |34| ) we can rewrite the general 
formula ( |l9| ) for the superfluid transition temperature in 
the form 



Tx = A—n 2 / 3 —F (b— 
m 770 V Vo 



rjo = fin 



where the function F(x) is defined by Eq. (J2C|) and 



A = 26, 



B = -ir 2 b£. 
3 



(35) 



(36) 



We see now that the critical temperature in this ap- 
proach is determined by the shear viscosity 77, and T\ 
will be smaller than the ideal gas value if the viscosity is 
considerably larger than 770 = fin, which may be called 
a 'quantum' of viscosity. For typical densities of liquid 
helium one has 770 — 2 • 10~ 5 poise, while the viscosity 
above the superfluid transition is 77 ~ 3.5 • 10~ 5 poise. 
Thus in the case of helium 77 is only slightly larger than 
770, and hence the temperature of the A transition should 
be close to that of the ideal gas, as it is experimentally. In 
sufficiently viscous Bose liquids the critical temperature 
of the possible superfluid transition should be suppressed 
and they will crystallize on cooling. 

The viscosity in helium slightly above T\ is almost 
independent of temperature and increases with density 
approximately by a factor of two when the pressure is 
raised from s.v.p. to 24 atm. We take experimental data 
for the pressure dependent viscosity at T = 3.2 K from 
Ref. |33| and turn it into n{n) dependence using the known 
equation of state at the same temperature. The density 
dependence of the ratio 77/770 is approximated then by the 
second order polynomial 



v/Vo-^Qin*) 1 



(37) 



!=0 



where n* — nr^ (ro — 2.556 A) is the reduced density 
and the coefficients are cq ~ 2.965, c\ ~ —15.701 and 
C2 — 34.027. Then we fit Eq. ( p5| ) to the experimental 
data for T\ in helium (actually to the Kierstead empirical 
equation for the A lineE-a) and obtain A ~ 1.39 and B ~ 
0.58, which correspond to quite reasonable values 



0.7, 



0.13 



(38) 



The resulting curve for the transition temperature is 
shown in Fig. 2 (solid line). We also extrapolate this 
curve to lower densities, corresponding to a metastable 
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liquid, though we do not expect our hydrodynamic for- 
mula ( [34|) to be valid at n ~ 0.2 -j- 0.3. At low densities 
the friction coefficient should be smaller than Eq. j3^ ) 
(see e.g. Eq. (|33]))and hence the true value of T\ should 
be higher. Also on this figure the condensation tem- 
perature for the ideal gas is shown (dashed line) along 
with the experimental data for liquid helium (squares) 
and recent PJMC results for metastable helium at nega- 
tive pressure!^! (circles). Note that the points at negative 
pressure were not used in the fit. 



R 2 ( T ) = (3/m/?)r(ft/3 - 



(40) 




0.20 
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0.30 



0.35 



0.40 



0.45 



0.50 



Density (nr 3 ) 



FIG. 2. Temperature of the superfluid transition versus 
the reduced density n" = nrf) (ro = 2.556 A). Solid line is 
the theory ( pH[ ) at A = 1.39, B = 0.58, squares denotes ex- 
perimental data for helium, dashed line corresponds to the 
ideal gap-. Recent simulation results for helium at negative 
pressureEj are also included (circles). 

The value of the Lindemann ratio £ from Eq. (|§8|) is 
quite close to that of the ideal gas £0 — 0.15 and also to 
£ ~ 0.12 obtained in Ref . [l3| by taking ^7 ~ K. If we now 
evaluate the value of the dampin g pa rameter 7 = A/m 
then at b ~ 0.7 we have from Eq. (34) fi-y ~ 13 K, which 
is indeed close to the zero point kinetic energy K ~ 15 
K. 

In the case of helium it appears possible to check also 
the validity of the Drude approximation (fTi] ) which re- 
sults in Eq. @ for the mean square displacement in 
imaginary time. For this purpose let us consider a more 
general quantity 

i? 2 (r) = ([r(r)-r(0)] 2 ) 

The R 2 discussed above is just the average of R 2 (r) over 
imaginary time. Proceeding as in Sec. II and taking the 
mobility in the Drude form ( |l7| ) we obtain 



OA + °° 

m M * J 



1 



./ > sin 2 (^) (39) 



At small r we can replace the sum in Eq. (^9|) by the 
integral to obtain R 2 (t) ~ (37L/m)r which is indepen- 
dent of 7 since on a very short time scale the particle 
moves without collisions. In the general case the sum 
can be calculated numerically and the result for (3 = 2 
K _1 and hj ~ 13 K (as obtained from the fitting of T\) 
is shown by the solid line in Fig. 3. Also the PIMC cal- 
culation results are shown as presented in Ref. [Thl where 
the same quantity was evaluated for bosons (triangles) 
and for distinguishable particles (squares). 




Note that for a free particle one would have 



FIG. 3. The diffusion of 4 He atoms as a function of imag- 
inary time r for T = 0.5 K(D = Ti/m and (3 = 1/T). The 
solid line results from Eqj— ^39[) at h-y — 13K. For comparison, 
PIMC simulation results!-] are also included, both for bosons 
(triangles) and distinguishable particles (squares). 

Is is clear from Fig. 3 that the Drude approxima- 
tion well account for the depression of R 2 (t) due to the 
interaction in the nondegenerate system. Since by our 
definition R 2 should be calculated for the system with 
all exchange effects neglected (see Sec. II) we conclude 
that the Drude approximation (and hence the Caldeira- 
Leggett model) may be quite sufficient for our superflu- 
idity criterion. 



V. 3 HE- 4 HE MIXTURES 

To check the universality of the Lindemann-like cri- 
terion it is necessary to consider some other superfluid 
systems. A natural choice to start with is the case of 
3 He- 4 He mixtures. One may try to repeat calculations 
from the previous section and express T\ in terms of vis- 
cosity. Unfortunately there seems to be no such viscosity 
measurements for a wide range of pressures as that of 
GoodwirH £m pure 4 He. Though there exists a recent 
investigatior£3 of rj in mixtures it deals only with pres- 
sures near the saturated vapor pressure. Above T\ (at 
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T ~ 3 K) the viscosity at s.v.p. is almost independent 
of temperature and gets smaller if the concentration x of 
3 He atoms increases. 

At small x this effect may be attributed simply to the 
decrease of the total density. Indeed, adding 3 He atoms 
to liquid 4 He leads to some reduction of the mean density 
(due to lighter mass of 3 He atoms they have larger zero- 
point energy and are harder to localize) and hence to a 
lower viscosity. But at a fixed density the local environ- 
ment of a given 4 He atom remains practically the same 
as for x — since interactions between atoms do not de- 
pend on the isotope mass. Hence the viscosity and the 
mean square displacement R 2 in dilute mixtures (small 
x) are likely to be the same as in pure 4 He at the same 
density. 

However, it is the mean distance between Bose parti- 
cles that enters in the r.h.s. of the criterion (w. This 
distance is now larger by the factor (1 — x)~ 1 <3 due to 
the presence of additional 3 He atoms which do not par- 
ticipate in superfluidity, so the criterion should read 

R 2 = ai-xy Va a 2 (41) 

where a = n~ 1//3 is determined by the total mean density. 

Now assuming that at a fixed density R 2 is indeed the 
same as in pure 4 He, and proceeding as in the previous 
section we obtain 

Tx = A^-n^F f—J—B^L) , m = hn 
m 770 \(l-x) z ' d r) J 

(42) 

where 77/770 is again given by Eq. (p7|). The universality 
hypothesis means that £ is independent of x, so that we 
may consider A and B to be defined by Eqs. (|3^) and 
( |38| ) . Then the only dependence of the superfluid tran- 
sition temperature on x is through the explicit factor of 
(1 - x)- 2/3 in Eq. @. 

To check the validity of the assumptions made we are 
to compare the formula (^2|) with, experimental data. 
Here we make use of experimentsLH where T\ was mea- 
sured in 3 He- He mixtures under pressure at several val- 
ues of x. We transform these results into T\{n) depen- 
dence using equations of state for different x at T = 1.5 
K, which hardly can introduce any significant error. 
These data are shown by squares in Fig. 4 for x = 0.089 
and x — 0.228 along with the pure 4 He results. Solid 
curves represent T\ as obtained from Eq. (|42). Note 
that only the x = curve was fitted to experiment. 

We see that qualitatively the formula ([|2]) correctly 
describes the dependence of the transition temperature 
both on density and on x, though for larger x's the the- 
oretical curves tend to lie lower than the experimental 
points. It is difficult to decide unambiguously whether 
this deviation indicate possible weak dependence of £ on 
x or simply the approximations made in evaluating R 2 
are too crude. In any case the estimate ( fi^ ) makes sense 
only for small x. 



2.4 



2.2 



1.6 



1.4 



- 




- 






x=0 




□ 






x=8.9% 












□ 






□ 








□ ^ 




x=22.8%^^^\ 








□ 

□ 






□ 






\ □ 









0.30 0.34 0.38 0.42 0.46 

Density (nr 3 ) 



FIG. 4. Critical temperature of 3 He- 4 He mixtures versus 
the reduced density n* = nrj) (ro = 2.556 A) for 3 He con- 
centrations x = 0, 8.9% and 22.8%. Solid line is the theory 
( [fil ) with A and B the same as for pure 4 He, while squares 
represent experimental data. Only the x = curve was fitted 
to experiment. 

VI. CONCLUSIONS 

In this paper the Lindemann-like superfluidity crite- 
rion proposed earlienl 3 ] was applied to several systems 
where the superfluid transition is known to occur, namely 
to a weakly non-ideal Bose gas, liquid helium and 3 He- 
4 He mixtures. The mean square displacement R 2 of a 
particle in imaginary time which enters in the criterion, 
and hence T\ were related here to the zero frequency 
friction coefficient (or to the shear viscosity in case of 
liquids). We found that in all cases the results obtained 
do qualitatively reproduce the observed behaviour of the 
transition temperature T\, as obtained from experiment 
(including numerical modeling). 

As far as the quantitative accuracy of the approach 
is concerned the situation is still not fully understood. 
It is clear, however, that one can hardly expect it to 
be very accurate. Indeed, the criterion proposed may 
be of real interest only if the critical temperature differs 
significantly from the ideal gas value (in the opposite case 
the estimate (|l|) is quite sufficient). But it is obvious e.g. 
from Eq. (|l6|), valid at Ti'y/To 3> 1, that the result in this 
case is rather sensitive to the exact values of the damping 
parameter 7 and the phenomenological parameter £ since 
they both stand in the exponential. These parameters are 
in general not known exactly and even a small variation 
in 7 or £ may significantly influence the resulting estimate 
forT A . 

It seems therefore that the best use of the formulas 
obtained may be the analysis of how the transition tem- 
perature, which e.g. is known for a given system, may 
change if some parameters of the system are slightly mod- 
ified. In this case one can fix all adjustable parameters 
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by f itting to the known data and then analyze their de- 
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viat ons for a modified system. A simple example of how 
this can be done was given in Sec. V where 3 He- 4 He 
mixtures were discussed. We have first fixed all free pa- 
rameters comparing the theory to the pure 4 He data and 
then take account of the nonzero concentration of 3 He 
atoms. 

The criterion proposed leads also to an unexpected pre- 
diction for the weakly interacting Bose gas. In this case 
the initial increase of T\ with interaction strength may 
be attributed to the change of the Lindemann parameter 
£ and this results in an upper bound on the peak value of 
T\/T®. Indeed, in the limit of zero density £ ~ 0.15 and 
at helium densities we have an estimate £ ~ 0.12 0. 13 
(see Ref. [l3] and Sec. IV) so that the change in £ and 
hence the maximum value of T\/T® could not be large. 
In Sec. Ill we obtain a rough .estimate T\/T® < 1.2 con- 
sistent with the PIMC resultaJ. Possible numerical eval- 
uation of £ in a dense system can make this inequality 
more accurate. 

Numerical estimates of £ would be desirable also from 
the general point of view. This is practically the only 
way to make a final decision concerning the validity of 
the approach presented here. One has to calculate R 2 
along the A line and then check whether £ = R 2 /a 2 can 
be regarded as a constant (at least in some density range) 
and to what extent is £ universal in systems with e.g. 
different interaction potentials. 
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usefull discussions and to D.M. Cepcrley for stimulat- 
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Fundamental Research (Project No 00-15-96698). 
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